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Polymer chain scission at constant tension - 

an example of force-induced collective behaviour 
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Abstract -The breakage of a polymer chain of segments, coupled by anharmonic bonds with 
applied constant external tensile force is studied by means of Molecular Dynamics simulation. We 
show that the mean life time of the chain becomes progressively independent of the number of 
bonds as the pulling force grows. The latter affects also the rupture rates of individual bonds 
along the polymer backbone manifesting the essential role of inertial effects in the fragmentation 
process. The role of local defects, temperature and friction in the scission kinetics is also examined. 
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Introduction. — The understanding of a great vari- 
ety of phenomena related to stability, fracture, and elastic 
behavior of materials requires fundamental knowledge of 
the inter molecular dynamics of bond breakage. In most 
cases scission of bonds may be caused by mechanical load, 
irradiation, or just increase in temperature. Examples re- 
lated to this field are diverse and include mechanical frac- 
ture of materials [T][5] , polymer rupture [2H7] , adhesion [8] , 
friction [2] , mechanochemistry [TUlfTT] . Recently, there has 
been an enormous increase of interest in polymer fragmen- 
tation due to the possibility of biomolcculc's micromanip- 
ulation in experiments using force spectroscopy methods 
(T2HH]-This has motivated also theoretical investigations 
and computer experiments [T5] . 

In particular, the problem of polymer fragmentation has 
got a longstanding history in scientific literature. The 
treatment of bond rupture as a kinetic process dates back 
to the publications of Bucche [TB] and Zhurkov et al. [T7J . 
In the recent years these seminal papers have been comple- 
mented by a variety of computer experiments. Molecular 
Dynamics (MD) simulations of chain rupture at constant 
stretching strain has been carried out, whereby harmonic 
[H1Q1], Morse [2DHH or Lennard- Jones [Hlj interac- 
tions have been employed. A theoretical interpretation 
of MD results, based on an effectively one-particle model 
(Kramers rate theory) has been suggested [55,26 . On the 
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other hand, an analytical treatment of a polymer fragmen- 
tation under constant stress have been proposed in terms 
of many-particle version of transition state theory |22j . 

Recently, we proposed a description of linear polymer 
scission under constant tensile force by using the multi- 
dimensional Langer-Kramers theory, which was found to 
compare favorably with the results of MD simulations [57] . 
Within this approach a single bond rupture is seen as a 
thermally activated escape from the bottom of a poten- 
tial well. The life time r before a bond scission takes 
place, is determined by diffusive crossing of an energy bar- 
rier Ef, that is reduced under the applied external force /. 
The adopted theoretical treatment assumes a single collec- 
tive unstable modes as being mainly responsible for chain 
breakage. Such unstable mode peaks around an "endan- 
gered" bond of negative spring constant and decays ex- 
ponentially towards both chains end. Similar collectivity 
effect has also been reported in the case of ring polymers 
stretched with constant strain |28j . 

In this letter we report some new results pertaining to 
the rupture kinetics of single ID and 3D polymer chains 
induced by constant tensile stress in a broad interval of 
pulling forces. 

The model. — As in our previous work [57], we use 
a coarse-grained model of a polymer chain of N beads 
connected by bonds, whereby each bond of length b is 
described by a Morse potential V M (r) = D{1 — cxp[— a(r — 



p-1 



J. Paturej et al. 



b)}} 2 , with a being a constant, a = 1, that determines bond 
elasticity, 

The dissociation energy D of a given bond is measured 
in units of fc^T, where ks is the Boltzmann constant and 
T denotes temperature. Since V* M (0) « 2, the Morse po- 
tential is only weakly repulsive and segments could par- 
tially penetrate one another at r < b. Therefore, in or- 
der to allow properly for the excluded volume interactions 
between bonded particles, we take the bond potential as 
a sum of V m and the so called Weeks-Chandler-Anderson 
(WCA) potential V WCA = 4e [(f) 12 - (f ) 6 + ±] 9(2 l / 6 a- 
r), with Q(x) = or 1 for x < or x > 0, and e = 1, a = 1. 
The parameter a sets up the length scale of equilibrium 
monomer size b = 2 1 / 6 er rj 1.12. The nonbonded inter- 
actions between monomers are also taken into account by 
means of the WCA potential. 

The dynamics of the chain is obtain by solving a 
Langevin equation for the position q n = [x n , y n , z n ] of 



each bead in the chain, mq„ = Fj! 



-7q„ + R„(i)- 



tnfinN (n, ■ ■ ■ , N), which describes the Brownian motion 
of a set of bonded particles whereby the last of them is 
subjected to external stretching force f = [/, 0,0]. The 
influence of solvent is split into slowly evolving viscous 
force and rapidly fluctuating stochastic force. The ran- 
dom, Gaussian force R„ is related to friction coefficient 7 
by the fluctuation-dissipation theorem. It should be noted 
that we consider only grafted chains, i.e. one of the chain 
ends is fixed in space. The integration step is 0.002 time 
units (t.u.) and time in measured in units of ^Jm/a 2 D, 
where m denotes the mass of the beads, m = 1. 

We start the simulation with all beads placed at distance 
b from each other, and then we let the chain to equilibrate 
in the Langevin heat bath. Due to the presence of the 
external pulling force, the equilibrium configuration of the 
chain is more or less stretched and deviates markedly from 
coil shape. Once equilibration is achieved, time is set to 
zero and one measures the elapsed time r before any of 
the bonds exceeds certain extension r^, which sets the 
criterion for considering such bond broken. We use a large 
value for the critical bond extension, = 56, which is 
defined as a threshold to a broken state. This convention is 
based on our checks that the probability for recombination 
of bonds, stretched beyond rh, is vanishingly small. We 
repeat this procedure for a large number of events 5 x 10 
so as to determine the mean rupture time (t) which we 
refer as Mean First Breakage Time (MFBT). The details 
of this method can be found in Ref . [27J . 

MD-results. — In our computer experiments we fo- 
cused on the following most salient properties of the bond 
breakage process: 

Chain length dependence of the MFBT (r) . In Fig.[lji) 
and the inset of Fig. [lji) we present numerical results for 
(t) as a function of the number of beads N for chains 
that are stretched in the interval 0.1 > / > 0.3 of pulling 
forces both in ID and 3D. Regardless of dimensionality 
of the examined systems, for a given value of / one ob- 
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Fig. 1: a) Mean first breakage time (r) vs. TV for a 3D chain. 
In the legend slopes of fitting lines (r) oc N^ 13 are presented 
which were found in the range N — 80-150. The inset shows 
(r) against N for 3D chains stretched by force / = 0.25 at 
different temperatures, b) Variation of slope /3 with external 
pulling force / for chains in ID and 3D. The inset shows 
(r) vs. N for a ID chain. Parameters of the heat bath are 
temperature T = 0.53 and friction 7 = 0.25. 



serves a power-law decrease, (r) oc N~P. This relationship 
is found for sufficiently long chains (asymptotic limit) - 
N > 80, where finite-size effects do not play a role. Addi- 
tionally, from the inset of Fig. QJi) we observe that there is 
no impact of the temperature on the value of the slope /3. 
Furthermore, Fig. [Tp) indicates that with growing tensile 
strength the life time (r) becomes nearly independent of N 
which is among the most important results of this study. 
This independence is fully consistent with our recent find- 
ings [27] for relatively strong pulling force, / = 0.25. It 
was shown in [37J that in this case the process of bond 
scission is governed by a collective unstable mode peaked 
around an "endangered" bond (i.e. a bond with negative 
spring constant) and decays exponentially towards both 
chain's ends. On the other hand, in the opposite limit of 
thermal degradation of polymers (/ = 0) (i.e., in the so- 
called thermolysis) the total probability for scission of a 
polymer with N bonds within a certain time interval is N 



p-2 



Scission of a polymer chain 



times larger than that for a single bond which is what one 
would expect if bonds do break entirely at random and 
independent of one another. The latter leads to the rela- 
tionship (t) oc TV -1 which has been seen recently in com- 
puter simulations of harmonic [25] and anharmonic |30j 
polymer chain models. Moreover, Fig. [Tp) clearly shows 
that with increasing pulling force / the exponent j3 grad- 
ually decreases within the interval < j3 < 1. Thus the 
slope /3 can be treated as a quantitative measure of the 
degree of coopcrativity in rupture events . As the slope f3 
decreases, the nature of scission events become more and 
more collective. 

It is pertinent to note that in the literature one finds 
conflicting data regarding the (r) vs. TV dependence. For 
example, in [26] it is claimed that (r) oc TV -1 , but a more 
close inspection of the Fig. 4 in this work where the depen- 
dence is shown gives for the exponent j3 = 0.17 -j- 0.2. In 
another work [28| . a chain rupture under fixed strain has 
been investigated by making use the many-body Langer- 
Kramcrs theory and MD-simulation. On the one hand, the 
authors discuss the presence of collective unstable mode 
(cf. [27] ) but on the other hand they claim that for the 
total chain again (r) oc TV -1 which should be considered 
as a hallmark of missing collectivity. Unfortunately, there 
have been no corresponding simulations which would have 
confirmed the conclusion made in 1281. 



Dependence of MFBT on Pulling Force. The depen- 
dence of the MFBT (r) on external force / for ID and 3D 
chains composed of TV = 30 beads is shown in Fig. [2] Ev- 
idently for sufficiently strong stretching forces / > 0.175 
an exponential decay (r) oc e^ Ea ~ a ^/ kBT is observed. The 
main reason of this is the following: As the pulling force 
grows, the energy barrier, which separates intact bonds 
from the broken ones, declines. As a consequence, (t) 
decreases. One should note that the parameters a and 
Eq change only slightly with the coupling parameter 7 of 
the thermostat. Fig. [2] indicates also a considerable dif- 
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Fig. 2: Force-dependent mean first breakage time for a ID and 
3D chains with TV = 30. The inset shows of (r) vs. Eb/T for 
ID system with T = 0.053 and 7 = 0.25. 

ference in the values of a between ID and 3D. In the 



inset of Fig. [2] we present (r) as a function of the ratio 
Ef, /T of the barrier height to temperature. This finding 
is in agreement with the understanding of the polymer 
rupture as a thermally activated process [161117) and is 
manifested by an Arrhenian relationship - (r) oc e Eb ^ T ', 
where E^ = Eq — af. 

Life-time probability distribution W{t). In Fig. [3] we 
display the probability distribution function W{t) of the 
observed scission times t for several ratios E^/T of the 




Fig. 3: Life-time probability distributions W(i) for different 
height of the energy barrier E^/T in ID. Here the chain length 
is TV = 30, the pulling force / = 0.15, and 7 = 0.25. Symbols 
denote simulation results and full lines stand for fitting func- 
tions W(t) oc exp(— t/r). The two insets show the dependence 
of t -1 on Eb/T [left panel], and on 7 [right panel]. 

barrier height to temperature in the case of ID chain 
composed of TV = 30 beads. It appears that W(t) goes 
asymptotically as W{t) oc e~ t l T ' yEh l T ^ in accordance with 
our recent findings [27]. Additionally from the inset [left 
panel] of Fig.[3]one may easily verify that the characteristic 
time goes as t(D{,/T) oc ex-p(EtfT). Moreover, a simple 
linear relationship has been found for friction dependence 
of t(j) as shown in the inset [right panel] of Fig. [3] 

Rupture probability histograms. In Fig. 0] we show the 
MFBT (t„) of the individual bonds for chains composed 
of TV = 30 beads in ID and 3D. We compare the results 
for different pulling forces as indicated in the legend. For 
the case / = 0.3 one may readily verify that the bonds lo- 
cated in the vicinity of the grafted bead live nearly twice 
as long as those close to the loose end where the tensile 
force is applied, regardless of dimensionality. In contrast, 
chains stretched by a gentle force, / = 0.15, display very 
different distributions of (t„). In the case of a 3D system 
the bonds at the tethered end live on the average signif- 
icantly (2-^3 times) shorter than those close to pulled 
end. Evidently, staring from n > 10 the lifetime of bonds 
progressively increases with increasing proximity to the 
free chain end. In the case of ID string stretched with 
/ = 0.15 this effect is missing and (r„) is uniformly dis- 
tributed along the polymer backbone apart from the first 



p-3 



J. Paturej et al. 



3D - f = 0.31 



A 

c 

V 



A 

c 

V 




5 10 15 20 25 30 

Consecutive Bond Number 



Consecutive Bond Number 



Fig. 4: Variation of the mean life time (r„) with consecutive 
bond number for a 30-particle chain in one and in three dimen- 
sions. The polymer is stretched by a force / = 0.15 or / = 0.3 
at T = 0.053 and 7 = 0.25. 



terminal bond which lives « 20% longer. As expected, all 
distributions displayed in Fig. [4] are asymmetric due to the 
constraint imposed on the motion of the first bead in the 
chain. 

In Fig. [5] we present the probability for bond scission 
of individual bonds in ID and 3D for several strengths 
of the pulling force. The histograms display the (normal- 
ized) rate at which a certain bond n along the polymer 
backbone breaks. From the inspection of Fig. [5^.) one sees 
that the preferential scission of the bonds with particu- 
lar consecutive bond number essentially depends on the 
value of force. For strong pulling / = 0.3 one finds that 
the terminal bond which is subjected to pulling as well 
as the bonds in its neighborhood break more frequently 
than whose around fixed end. Evidently, in this case the 
rupture rate decreases steadily from the free chain end to 
the tethered one. A similar scission scenario is visible also 
for the ID chain as shown in the inset of Fig. [5^,). 

In contrast, as the stretching force is decreased, the 
corresponding rupture histogram for a 3D chain becomes 
flatter. For / = 0.2 the distribution of scission rates be- 
comes uniform exept for the bonds in the vicinity of both 
ends. A further decrease of the pulling force results in 
a qualitative change in the distribution. Evidently, for 
/ = 0.15 the bonds in the middle of the chain, which are 
also somewhat closer to the fixed chain end, become more 
vulnerable as compared to those at the chain ends. Note 
that for the smallest pulling force (/ < 0.15) the rupture 
histogram already resembles the respective histogram in 
the case of thermal degradation of a polymer [50] which 
takes place in the absence of externally induced tension. 
If two pulling forces are applied simultaneously to both 
chain ends, Fig. [5p)-d), one finds expectedly a symmetric 
scission probability distribution regarding bond number, 
even though (r) hardly changes - Fig.U). Again, the scis- 
sion rate of the terminal bonds goes down with decreasing 
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Fig. 5: a) Rupture probability histograms for ID and 3D chains 
composed of N = 30 for different pulling forces as indicated. 
b)-d) Rupture probability histograms for 3D chains made of 30 
particles pulled at both ends with a tensile force /. e) Variation 
of (t) with force / for chains pulled either at the one end, or 
at both ends. Here T = 0.053 and 7 = 0.25. 



strength of the pulling force /, and for / = 0.15 a local 
maximum in the rate builds up for the bonds that are in 
the middle of the chain. 

How can such an inhomogeneity in the probability of 
bond rupture be understood? A possible explanation of 
the change in the location of preferential breakdown sites 
along the chain may be gained by Fig. [5] In this figure we 
present maps of the density distribution P(x,r) of bead 
positions where x is measured in direction of the pulling 
force / whereas r = y y 2 + z 2 denotes the radial compo- 
nent. Fig. [5] [upper panel] indicates that at high stretch- 
ing (/ = 0.3) the most probable position of the beads is 
along the direction of the tensile force. For / = 0.3 (strong 
stretching) the chain conformation corresponds to a quasi- 
1D structure, and the transversal fluctuations are reduced. 
In contrast, when the pulling force is weak (/ = 0.15) 
one finds from Fig. |5J [lower panel] that the individual 
beads are free to make big excursions in space - P(x,r) 
is roughly two times broader in the middle. Thus, Fig. [S] 
suggests that the density maps comply with the rupture 
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Fig. 6: Probability density distribution P(x, r) of beads in a 3D 
chain with N = 30 particles at force: / = 0.3 [upper panel], 
and / = 0.15 [lower panel]. The x-axis coincides with the 
direction of pulling force whereas r = ^/y 2 + z 2 denotes radial 
component of the bead position. Different colors indicate the 
value of the PDF as indicated in the legend. Here T = 0.053. 



histograms given in Fig. [5£i). For / = 0.3, due to larger 
freedom around the pulled end, the end bonds break more 
easily in the terminal part of the chain. When the force 
is weak, / = 0.15, the beads become more mobile around 
the center of the polymer which in turn leads to increased 
bonds scission rate there. 

Chain Defects. In Fig. [7] we present the results of 
computer experiments concerning rupture of 30-particle 
chains in which a single defect is introduced. We fo- 
cus on two kind of defects. First we examine the effect 
of an "isotope" -like defect in which a mass rrid of the 
monomer located in the center of a chain is changed, while 
the masses of all remaining beads remain unchanged and 
equal to m = 1. In Fig. [7jt) we compare the rupture 
histogram for a chains with such a heavy /light bead to 
the respective scission probability distribution in a uni- 
form system. If the central bead is replaced by a parti- 
cle which is lighter /heavier than the rest of the segments, 
the rupture probability increases/decreases in the imme- 
diate vicinity of the introduced mass-defect. This effects a 
small group of beads ~ 3 and is pronounced only when the 
stretching force is sufficiently small (/ < 0.2). Evidently, 
a lighter particle can be kicked more easily by the thermo- 
stat whereby inertial effects would stretch both adjacent 
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Fig. 7: a) Rupture probability histograms for 30-bead chains 
with a single defect bead of mass mj introduced in the center 
of a chain. Here / = 0.15. The inset shows the mean life time 
time dependence on the mass of the defect, b) Rupture prob- 
ability histogram for chains with a single bond defect (bond 
strength Da) in the middle of the chain. The inset shows de- 
pendence of mean first breakage time on the strength of the 
defected bond. Parameters of the heat bath are: T = 0.05328 
and 7 = 0.25. 



bonds beyond the scission threshold. In contrast, as is seen 
from the inset of Fig. UK), the MFBT (r) grows as the de- 
fected bead becomes heavier. A defect with a larger mass 
is hard to accelerate, its mobility is low, and its bonds 
remain unstrctched. In some interval of time such de- 
fect experiences many kicks which effectively cancel each 
other before a dangerous bond stretching occurs. Thus, 
the chain becomes locally more immune against breakage 
events which is reflected by the decrease in the probability 
of rupture. 

As expected, the tensile strength of a particular bond 
(measured in units of D) affects the scission process too. 
To see this we varied the bond strength Dd of the middle 
bond. In Fig.[7p) we present rupture histograms for chains 
with N = 30 and different values of Dd- Clearly, the small 
variation of bond strength (±10%) results in a change of 
the scission probability which is located exactly at the 



p-5 



J. Paturej et al. 



defect position. In the inset of Fig. [7p) we present the 
MFBT (r) function of Dd- One can see that the 

introduction of weaker bonds Dd < 1 in the chain results 
in a decrease of (r). This is due to the fact that the 
lifetime td for the weak bond is very short and therefore 
dominates the mean (r) . On the contrary, it appears that 
(t) is not sensitive for defects with Dd > 1. In this case 
defected bonds appear to resist scission events which then 
happen predominantly in the remaining bonds. 
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Concluding remarks. 

marized as follows: 



Our findings can be sum- 



• The mean life time of the polymer chain at constant 
tensile force depends on chain length like (r) oc N~p 
whereby the power law exponent /? varies in the inter- 
val < (3 < 1. Generally, it appears that the expo- 
nent /3 systematically declines as the external pulling 
force / grows. This behaviour indicates a growing 
degree of coopcrativity during the chain breakage as 
the pulling force / is increased. 

• The MFBT follows an Arhenian law (r) oc exp(E b /T) 
whereby the effective activation barrier for scission 
changes with varying pulling force / as Ef, = Eq — af 
in line with earlier theoretical predictions [T7]. The 
scission times in a polymer chain under tension are 
exponentially distributed, W(t) oc exp(—t/T(Eb/T)). 

• The rates of bond rupture are distributed differently 
along the polymer backbone in the ID and 3D chain 
models. In a ID chain the rupture rate steadily grows 
as one approaches the free chain end where the ex- 
ternal pulling force is applied whereas in a 3D chain 
bonds break predominantly in the middle of the chain. 
Bond rupture histograms correlate with the degree of 
spreading in the monomer density distribution, indi- 
cating that scissions occur most frequently in those 
parts of the macromolcculc which undergo large fluc- 
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ity provide a plausible interpretation of the observed 
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• The probability histograms for bond scission provide 
a clear picture of the impact of polymer defects on 
the fragmentation process and underline thereby the 
role of inertial effects. Bonds, connected to lighter 
segments break become a preferred site of rupture 
whereas heavier segment stabilize the chain. 
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